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Abstract 


We have studied the contribution of the high-scale SUSY to the Kl —y and K + —> 

processes by correlating with the CP violating parameter ex- Taking account of the 
recent LHC results for the Higgs discovery and the SUSY searches, we consider the high- 
scale SUSY at the 10 — 50TeV scale in the framework of the non-minimal squark (slepton) 
flavor mixing. The Z penguin mediated the chargino dominates the SUSY contribution 
for these decays. At the lOTeV scale of the SUSY, the chargino contribution can enhance 
the branching ratio of Kl —>• in eight times compared with the SM predictions 

whereas the predicted branching ratio BR{K + —>• tt + uv) increases up to three times of 
the SM one. The gluino box diagram dominates the SUSY contribution of ex up to 
30%. If the down-squark mixing is neglected compared with the up-squark mixing, the Z 
penguin mediated the chargino dominates both SUSY contributions of BR(Kl —>• iT°nu) 
and ex- Then, it is found a correlation between them, but the chargino contribution to 
ex is at most 3%. Even if the SUSY scale is 50TeV, the chargino process still enhances 
the branching ratio of Kl from the SM prediction in the factor two, and ex is 

deviated from the SM prediction in 0(10%). We also discuss the chargino contribution 
to Kl 7r°e + e _ process. 
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1 Introduction 


The K meson physics have provided important informations in the indirect search for New 
Physics (NP). Especially, the rare decay processes K + — y tt + uu and Kl —> are known 

as the clean one theoretically DEI- Therefore, these both processes have been considered to 
be one of the powerful probes of NP El-1131 whereas these decay widths are bounded by so 
called the Grossman-Nir bound for the NP [T5l fT6] . 

The Kl —y process is the CP violating one and provides the direct measurement 

of the CP violating phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix [iT7l [IB]. In 
addition, the CP conserving process K + —y tt + uu is also the physical quantity related with 
the unitarity triangle (UT). On the other hand, the CP violating parameter ek, which is 
induced by the K° — K° mixing, also constrains the height of the UT. Hence these measured 
variables give us the information of the UT fit as well as the CP violating quantity sin 2(f) \ 
induced by the B° — B° mixing. Furthermore, the K —y nvv processes are expected to open 
the NP window in the CP violating flavor structure. 

The K + — y n + vv and Kl — y ir°vv decay processes are governed by the Z penguin diagram 
in the Standard Model (SM) [T9], which predicts 

BR(K l ->■ vr 0 z/h) SM = (2.43±g;|? ± 0.06) x 10" 11 , 

BR(K+ -)• n + uP) SM = (7.81±g;f? ± 0.29) x 10' 11 . (1) 

In the estimation of the branching ratio of K —>■ nvu, the hadronic matrix elements can be 
extracted with the isospin symmetry relation [20, [2Tj . These processes are theoretically clean 
because the long-distance contributions are small [T2] . and then the theoretical uncertainty 
is estimated below several percent. On the other hand, ek has the different flavor mixing 
structure from these processes since it is induced by the box diagram of K° — K° mixing. 
Therefore, the NP is expected to appear in both K —> nvv and Ek with different magnitudes. 

On the experimental side, the upper bound of the branching ratio of Kl —> 7T°z/h is given 
by the KEK E391a experiment [22]. The branching ratio of K + —> n + vD measured by the 
BNL E787 and E949 experiments is consistent with the SM prediction [S3] ; 

BR(K l vr 0 i/p) ex p < 2.6 x 10~ 8 (90%C.L.), 

BR(K + ->■ 7 T + vv) exp = (1.73i};^) x lO' 10 . (2) 

At present, the J-PARC KOTO experiment is an in-flight measurement of Kl 7i 0 uu 
approaching to the SM predicted precision [24., [22], while the CERN NA62 experiment [26j 
studies the K + —y n+vi' process. 

On the theoretical side, the supersymmetry (SUSY) is one of the most attractive candi¬ 
dates for the NP. However, the SUSY signals have not been observed yet, and then the recent 
searches for new particles at the LHC give us important constraints for the SUSY. Since the 
lower bounds of masses of the SUSY particles increase gradually, the squark and the gluino 
masses are supposed to be at the higher scale than 1 TeV [2711281 [29] . Moreover, the SUSY 
models have been seriously constrained by the Higgs discovery, in which the Higgs mass is 
126 GeV [301I3JJ. 

These facts suggest a class of SUSY models with heavy sfermions. If the squark and 
slepton masses are expected to be 0(10 — 100) TeV, the lightest Higgs mass can be pushed 
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up to 126 GeV, whereas all SUSY particles can be out of the reach of the LHC experiment. 
Therefore, the indirect search of the SUSY particles becomes important in the low energy 
flavor physics [32] [33, 33] . 

So far, the effects of SUSY on the K + —)■ n + vv and Kl —)■ processes have been 

studied in the framework of the Minimal Supersymmetric Standard Model (MSSM) with the 
minimal flavor violation (MFV) scenario intensively [8, TO]. Since the SUSY mass scale is 
pushed up more than 1 TeV region at present, the effect of the MSSM with MFV is expected 
to be very small. These processes are also discusses in the framework of the general SUSY 
model pj '33 021 ST, ESI ESI 3D] at the (9(500) GeV scale. 

We have studied the SUSY contribution to the CP violation of the B meson and ex 
induced by the K° — K° mixing under the relevant SUSY particle spectrum constrained by 
the observed Higgs mass [13]. Then, it is found that the SUSY contribution could be up to 
40% in the observed ex, on the other hand, it is minor in the CP violation of the B meson at 
the high scale of 10 — 50 TeV. Therefore, in this paper, we investigate the high-scale SUSY 
contribution to K + —> n + vv and Kl —)■ n Q vv by correlating with ex in the framework of the 
mass eigenstate of the SUSY particles, which is consistent with the updated experimental 
situations like the direct SUSY searches and the Higgs discovery, with the non-minimal squark 
(slepton) flavor mixing. 

Our paper is organized as follows. Sec.2 gives the basic framework of K + 7r + uu, 
K l —> and ex in the SM and the MSSM. In Sec.3, we present the setup of the high- 
scale SUSY. In Sec.4, we discuss our numerical results. Sec.4 is devoted to the summary. 
The SUSY mass spectra and the Z penguin amplitude mediated the chargino are given in 
Appendices A and B, respectively. 

2 Basic framework 

In this section, we present the basic formulae for the K —)■ nuD decay and the CP violating 
parameter ex, which correspond to |A S\ = 1 and |A S\ = 2 processes, respectively. The 
K + —>■ n + vv and Kl —>■ n°vv processes are clean ones theoretically since the hadronic matrix 
elements can be extracted including isospin breaking corrections by taking the ratio to the 
leading semileptonic decay of K + —>■ 7i 0 e + is. Moreover, the long-distance contributions to 
these rare decays are negligibly small. Therefore, the accurate measurements of these decay 
processes provide the crucial tests of the SM. Especially, the Kl —> process is purely 

the CP violating one, which can reveal the source of the CP violating phase. 

On the other hand, the CP violating parameter ex is measured with enough accuracy. 
The major theoretical ambiguity comes from the hadronic matrix element factor B x- The 
recent lattice calculations give us the reliable value for Bx m 03 . The more accurate 
estimate of the SM contribution enables us to search the NP such a SUSY because we know 
the accurate observed value of ex- Actually, the non-negligible SUSY contribution has been 
expected in ex at the scale of (9(100) TeV [02] [331 S3] • Consequently, it is required to examine 
the high-scale SUSY contribution in K —>■ nvv by correlating with ex- 
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2.1 Basic framework : K + —>■ and Kl -f iT°isi> 

2.1.1 K + —>■ 7T + z/i? and Kl —)• 7r°^P in the SM 

Let us start with discussing the framework of the A" + —» n + vv and Kl —> 7T°z/h processes in 
the SM [I], The effective Hamiltonian for K —> nvV in the SM is given: 

%r E KUiVc + V£VWf.] (UtX) + H.c., (3) 

V2 7TS111 0 W . ^ 
v i=e,/i,r 

which is induced by the box and the Z penguin mediated the W boson. The dominant 
box contrition is derived by the top-quark exchange, on the other hand, the charm-quark 
exchange contributes to the Z penguin process as well as the top-quark one. The up-quark 
contribution is negligible due to its small mass. So, the loop function X c denotes the charm- 
quark contribution of the Z penguin, and X t is the sum of the top-quark exchanges of the 
box diagram and the Z penguin in Eq. (lHl) . 

Let us define the function F as follows: 

F = V'V^X' + v;,v ti x, . (4) 

The branching ratio of K + —> 7i + uu is given in terms of F. Taking the ratio of it to the 
branching ratio of K + —> 7T°e + h, which is the tree level process, we obtain a simple form: 

BR(K + —> 7 t + vv) 2 / a 

BR(K+ 7i°e+u) = \V US \ 2 V27rsin 2 0vF 

The K + —> n°e + u decay is precisely measured as BR(K + —> 7r°e + z/) exp = (5.07±0.04) x 10“ 2 
[43] . and its hadronic matrix element is related to the one of K + —> 7r + vv with the isospin 
symmetry: 



(vr°| (d L ^s L ) \K°) = (tt 0 | {slI^ul) \K + ), (6) 

(tt + | (slY^l) | K + ) = v / 2(7t°| (s l Yul) \K + ), (7) 

where the coefficients are determined by the Clebsch-Gordan coefficient. By using this rela¬ 
tion, the hadronic matrix element has been removed in Eq.(l5l). 

Now the branching ratio for K + —> n + vD is expressed as follows: 

BR(K + ->■ tt + uP) = 3k • r K+ \F\ 2 , (8) 


K = 


IK. 


a 


27rsiir6 ) iy 


BR(K + ->• 7r°e + h), 


(9) 


where r^+ is the isospin breaking correction between K + —>■ 7T°e + P and K + —> 7T°e + P [20ll2l] . 
and the factor 3 comes from the sum of three neutrino flavors. It is noticed that the branching 
ratio for K + — >■ 7 t + uu depends on both the real and imaginary part of F. 
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For the Kj j —y n Q i>v decay, the K° — K° mixing should be taken account, and one obtains 
A(K l ir°vv) = 2a 2 {y l L ~iy L ) (tt 0 | [F(s L ^d L ) + F*(d L 7 M s L )] | K L ) 

v 2 7TS111 U\,y 

= 2(X 2 n ( p l7/Fl) ~y= [^(! + e)(n°\(s L lM\K 0 } + F*( 1 - e)(7r°|(J L 7 /i s L )|/l 0 )] 

V 2 vrsin o\y y 2 

= % 2< ?n -7= t F ( X + e) - F*(l - e)] (tt°| (dL^SL^K 0 ) 

V 2 7 TSin Oyy Y 2 

" (^X) 75 2ImF XI(Xv7>r>. < 10 ) 

In the step of the first line going to the second line in (TTUT) , we use 

I K l ) = F [(1 + e)\K°) + (1 - e)|A-»>] , (11) 

and then, after using the CP transition relation in the second line, 

CP|F°) = -\K°), C| K°) = \K°), (12) 


(n°\(d L ^s L )\K°) = -(n°\(s Ll ,d L )\K°), (13) 


we obtain the equation in the third line. In the final line, we neglect the CP violation in 
K 0 — K° mixing, e, due to its smallness |e| ~ 10~ 3 . 

Taking the ratio between the branching ratio of K + — y n°e + v and Kl —t 7r°KK, we have 
the simple form: 


BR(K L ^*°vu) 2 f a \ 2 t(K l ) v 2 
BR(K+ -> n°e+u) \V US \ 2 \2nsm 2 9 w J t(K+) r 


(14) 


Therefore, the branching ratio of Kl —>■ tt°uv is given as follows: 

BR(K l 7r°i / p) = 3 k • !^liM(ImF) 2 , (15) 

r K + r(A+) 

where tk l and tk+ denote the isospin breaking effect [ 20l l2Tj . It is remarked that the 
branching ratio of Kl —>■ depends on the imaginary part of F. Since the charm-quark 

contribution is negligible due to the small imaginary part of V* s V c d, it is enough to consider 
only the top-quark exchange in this decay. 

In the SM, K + —y n+vv and Kl —>• ^t°vis are related to the UT fit. We write down the 
branching ratio in terms of the Wolfenstein parameters. Since ReF and ImF are given as 


ReF = —XX C - A 2 A 5 (l - p)X t , ImF = A 2 X 5 r)X t , (16) 


we can express the branching ratio of these decays as 


BR(K + -»■ 7 x + uu) 


3K-r^+[(ReF) 2 + (ImF) 2 ] 


3k • r K+ ■ A 4 X 10 X 


(P 



(17) 
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where 


Y c 

P ° ~ + A 2 X a X, ' 


( 18 ) 


and 


BR{K l -> vr °uu) = 3k • TKl ' (ImF) 2 

r^+ r(A+) 

< 19 > 

It is noticed that BR(K + —>• 7r + z/i/) in Eg (ITT)) is approximately a circle centered at p = p 0 ~ 
1 .2, fj = 0 on the p-fj plane. On the other hand, BR(Kl —y tt 0 uu) in Eg ([T9]) just depends on 
r) and it can determine the height of the UT directly. In this way, the precise measurements 
of K + —> n + vv and Kl —> become crucial tests for the SM. 

Before going to discuss the SUSY formulation, we present the general bound between 
K + —> n + vD and K L —y i t°uD, so called the Grossman-Nir bound JT5]. As seen from above 
formulations, since the two processes are determined by the imaginary part and the absolute 
value of the same coupling, the model independent bound is obtained: 

BR(K l 7r°m7) < r * L ■ BR(K + ->• n + vv) < 4.4 x BR(K + ->• ir + vv) , (20) 

where we use the isospin symmetry A(K + —> n + uu) = \/2A(K 0 —> n°uP). This bound must 
be satisfied for any NP (15, 16] . 

2.1.2 K + —> ti + vv and Kl —> 7r 0 uP in the MSSM 

The effective Hamiltonian in Eg. (|3]) is modihed due to new box diagrams and penguin dia¬ 
grams induced by SUSY particles. Then, the effective Lagrangian is given as 

£eff= v [(4, (sA) + c*4, Wi)] tbl) + H.C. , (21) 

where i and j are the index of the flavor of the neutrino final state. Here, C*y LL VRL is the 
sum of the box contribution and the Z penguin one: 

rip — _ p 21 *j _ a YM p21 rij 
°VLL — -°VLL ZL° : 

rK _ _p 21 *i _ a<2 p 21 yo fpo) 

WRL — - d VRL ZR° ) 

where the weak neutral-current coupling Qzl = 1/2, and B^f^ h and P|^(r) denote the box 
contribution and the Z penguin contribution, respectively. The V, L and R denote the vector 
coupling, the left-handed one and the right-handed one, respectively. In addition to the W 
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boson contribution, there are the gluino g, the chargino x ± and the neutra lin o x° mediated 
ones We write each contribution as follows: 

Bvl J l = (x*) + <i(x°) , 

tx*)+<i(x°), 

Pit = P&W) + P&H±) + + Pztb?) + fS2(x") , 

P£k = fzlte) + -Pzttx*) + Pzr(x 0 ) , (23) 


where (i,j) denotes the neutrinos of final state. Explicit expressions are given in Ref. [44]. 
It is well known that the most dominant contribution comes from the Z penguin mediated 
chargino for the K + —y ix + vv and Kl —> tt°uu decays [12]. 

The branching ratio of K + — y n + i>v and Kl —)■ are obtained by replacing internal 

effect F in Eqs. (JSJ) and (JT3J) to Cy LL + Cy RL as follows: 


BR(K* -» x+ V u) = K ■ r If+ V |C« LL + CVI 2 

i=e,fi,r 


BR(K l -f n°uu) 


= K • 


r K L t(K L ) 
r K + t(K+) 


^ |Im(C'y LL + C'vrl)! 2 ■ 


(24) 

(25) 


2.2 ek in the MSSM 

It is well known that the CP violating parameter Ek induced by the K° — K° oscillation gives 
us one of the most serious constraint to the NP. The general expression for ek is given as 


€k 



f lm (Mg) 
V A M k 



ImA 0 
ReA 0 ’ 


(26) 


where A 0 is the 0-isospin amplitude in the K —>■ nn decay, and My 2 is the dispersive part 
of the K° — K° oscillations, and AM K is the mass difference of the neutral K meson. The 
effects of £ ^ 0 and < 7 t/4 were estimated by Buras and Guadagnoli [45] . In the SM, the 
off-diagonal mixing amplitude M(\ is obtained as 


M 


12 

K 


(K 0 \Has=2\K°) 

5 {^) 2 M wBkFI m k 


%c(V c ,V^) 2 S(x c ) + r,u(V„V£) 2 S(x t ) 


+ 2r ta (V c ,V c - d )(V t ,V^)S(x c ,x t ) , 


(27) 


where S(x) denotes the SM one-loop functions [46], and rj C c,tt,ct are the QCD corrections [45] . 
Recent lattice calculations give us the precise determination of the B K parameter mm- 
Once taking account of the NP effect, the expression of M[ 2 is modified. In the case of 
the SUSY, new contributions to the box diagrams are given by the gluino g, the charged 

1 The wino-higgsino mixing is tiny in our mass spectrum. 
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Higgs H ± , the chargino x* 1 and the neutralino x° exchanges: 


M* = M/f M + m£ susy 

= M*{W) + M*{H±) + M^(x ± ) + M«(x°) + M«{g) + M«( X °g). 


The explicit formula has been presented in Ref. [ 44] . 


3 Setup of the squark flavor mixing 


We present the setup of our calculation in the framework of the high-scale SUSY. Recent 
LHC results for the SUSY search may suggest the high-scale SUSY, 0(10 — 1000) TeV 
|32l |33l IM1 HU since the lower bounds of the gluino mass and squark masses exceed 1 TeV. 
Taking account of these recent results, we consider the possibility of the high-scale SUSY at 
10, 50 TeV, in which the K —> i tuu decays and are discussed.. 

Another important experimental result should be mentioned is the Higgs discovery. The 
Higgs mass mn — 126 GeV gives effect to the SUSY mass spectrum. In general, there are 
two possibility to get Higgs mass value, one is the heavy stop around 10 TeV, and the another 
is the large X t = A 0 — /icot/3 given by the A-term. In the case that the SUSY scale is 10 
to 50 TeV, we have already obtained the SUSY mass spectra which realize the Higgs mass 
at the electroweak scale with Renormalization Group Equation (RGE) running in previous 
work [34]. We use this numerical result for the SUSY particle mass spectrum. In this study, 
the 1st and 2nd squark are almost degenerated due to the assumption of the universal soft 
masses. On the other hand, the 3rd squark mass obtains the large contribution from the 
RGE running due to the large Yukawa coupling of the top-quark. Therefore, the mixing 
between 1st and 2nd is negligible, and it is taken account in the subsequent discussion for 
squark flavor mixing. The SUSY spectra at 10 and 50 TeV are given in Appendix A. 

Once the SUSY mass spectrum is fixed, we can calculate the left-right mixing angle 9 q , 
which is defined as 


d m b (A 0 — /rtan/3) 

9 ~ -5-5- 

m~ — m~ 

i>L b R 


d u ~ 


m t (A 0 — fi cot /?) 

9 2 

m~ — m~ 

+L tR 


(28) 


In the case of the SUSY scale to be 10 and 50TeV, the left-right mixing angles of squarks 
and sleptons are very small as ( 6 d ~ 0.0062, 6 U ~ 0.0024, 0 e ~ 0.014) and (6 d ~ 0.0009, 
9 U ~ 0.0007, 9 e ~ 0.005), respectively. 

The SUSY brings the new flavor mixing through the quark-squark-gaugino couplings and 
the lepton-slepton-gaugino ones. The 6x6 squark mass matrix M| in the super-CKM basis 
turns to the mass eigenstate basis by diagonalizing with rotation matrix T^ as 


m\ = T (9) M|r {<z)t , (29) 

where T^ is the 6x6 unitary matrix, and we decompose it into the 3x6 matrices as 
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rfo> = (T^\ T (q)) ) T in the following expressions: 



( 


c qL 

c 13 

0 

sf R e *^13 Cgq 

0 

0 

—sf^e 3 s eq e l P 

\ 

p(<?) 

1 L ~ 


s qL 

-’23 

s^e'C^i 3 -^ 23 ) 

c qL 

c 23 

stc^e-^cgq 

0 

0 - 

' s 23 c 13 e_ *^ 3 SgqP’P 



V 


s gcge< - 

■s^e^3 

qL qL 
c 13 C 23 C d q 

0 

0 

Scrips 



^0 

0 

s qR seie~ 1 ^ 3 e 

-i(p q 

r qR 

c 13 


0 

s qR e-< 3 c d q 

\ 

p(9) _ 

1 R ~ 

0 

0 





r qR 

c 23 

q q R r qR p -i<t>%3 rn 
s 23 C 13 e c 9 q 


\o 

0 

qRqR -i<yJ 

c 13 c 23 s 0 q & 

qR qR id> q ^ 
-«13 C 23 e 013 


qR iSld} qR qR 

Sh e - C 13 C 23 C 0 q 

/ 


(30) 

where we use abbreviations c qR,qR = cos 9 qR,qR , s q ^ qR = sin 8 qR ' qR , cg q = cos 9 q and sg q = 
sin# 9 . It is remarked that we take s q Rqh = 0 due to the degenerate squark masses of the 
1st and the 2nd families as noted in Appendix A. The angle 6 q is the left-right mixing angle 
between and qn, and they are calculable as mentioned above. Then, there are free mixing 
parameters 9 q ^ ,qR and <pjf’ qR - For simplicity, we assume s q j = s qR . On the other hand, we 
scatter (p qR and (j) qR in the 0 ~ 27T range independently. It should be noted that the mixing 

angles s q ^ R ^ have not been constrained by the experimental data of B, D and K mesons in 
the framework of the high-scale SUSY [33] . 

For the lepton sector, the mixing matrices T^^ have the same structure as the quark 

one in the charged-lepton flavors, however, there is only the left-handed in neutrinos. 

As well known, the charged Higgs and the chargino contributions dominate the K —y nvv 
processes [12]. Since the SUSY scale is high in our scheme, the charged Higgs are heavy, 
d(10TeV), so the charged Higgs contribution is suppressed in our framework. On the other 
hand, the dominant SUSY contribution to ex comes from the gluino box diagram if the 
flavor mixing angles of the down-squark and the up-squark are comparable. In addition, the 
chargino box diagram is also non-negligible. Consequently, we will discuss the both cases in 
which the down-squark mixing angles are negligible small and are comparable to the 

up-squark mixing angles . We scan the phases of Eq. (j3U|) for up-squarks, down-squarks, 
charged-sleptons and sneutrinos in the region of 0 ~ 27T independently. 

In our framework, the K —y nuu processes are dominated by the Z penguin mediated 
the chargino exchange, Pzt(x ±S ) i n Eq.([23]) , which are occurred through the tLS^d^X^ and 
interactions, respectively. In our basis, the relevant mixing is given by 

tUir = (ri”V CK M);(t f +)f+-(UYvfcKMWfv)?. pi) 

92 

where q = s,d, I = 1 — 6 for up-squarks, and a = 1,2 for charginos. The Vckm is 
the CKM matrix, and U+ is the 2 x 2 unitary matrix which diagonalize M' c Mc, where 
Me is the 2x2 chargino mass matrix. The fu denotes the yukawa coupling defined by 
frjv sin /3 = diag (m u ,m c ,m t ). Therefore, the combinations of mixing angles and phases in 
Eq. (j30]h c ?3 s ?3 S 23 e*^i 3 ) and c\ R s\ R sl R e 1 ^ 3 ) are important for our numerical analy¬ 
ses in the next section. We show the formula for P^tix^) i 11 Appendix B. 






-J 

c? 



BR(/f + -*7r + vv) 


Figure 1: The predicted BR(K L —» n°uD) versus BR(K + —» n + vv) at the SUSY scale of 10 
TeV with the mixing angle of s u = s d = 0.1. The pink cross denotes the SM predictions. The 
red dashed lines are the la experimental bounds for BR(K + —> n + vv). The green slanting 
line shows the Grossman-Nir bound. 

4 Numerical analysis 

Let us discuss the high-scale SUSY contribution to the K —y ttuD processes by correlating 
with 6k P3j- At present, we cannot confirm whether the SM prediction efU is in agreement 
with the experimental value e^ p because there remains the theoretical uncertainty with an 
order of a few ten percent. However, the theoretical uncertainties of ex are expected to 
be reduced significantly in the near future. Actually, the lattice calculations of B K will be 
improved significantly [4T1 142] . whereas \V c b\ and the CKM phase 7 will be measured more 
precisely in Bclle-II. Therefore, we will be able to test the correlation between K — y nuD and 
£k- 

In our previous work, we have examined the sensitivity of the high-scale SUSY with 10 
and 50 TeV to ex- It is found that the SUSY contribution to ex is allowed up to 40%. We 
begin to discuss the SUSY contribution at the 10 TeV scale. The present uncertainties in 
the SM prediction for ex are due to the CKM elements V c b, p and fj, and the Bx parameter. 
We take the CKM parameters V c b, p and fj at the 90 % C.L. of the experimental data: 

\Vcb\ = (41.1 ±1.3) x 10" 3 , p = 0.117 ±0.021, fj = 0.353 ± 0.013. (32) 

For the Bx parameter, the recent result of the lattice calculations is given as mn 

B K = 0.766 ± 0.010 , (33) 

which is used with the error-bar of 90% C.L. in our calculation. 

In the beginning, we show the numerical results at the SUSY scale of 10 TeV. FigJT] shows 
the predictions on the BR(Kl —» 7 T°j/h) vs. BR(K + — y 'K + uv) plane, where phase parameters 
are constrained by the observed \ex\ with the experimental error-bar of 90%C.L. Here, we fix 
the mixing parameters in Eq. (j30l) by taking the common value s^ = s u — 0.1 (i — 1, 2) 
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Figure 2: The predicted BR(K L —> n°vv) versus the SUSY contribution ratio of ex at the 
SUSY scale of 10 TeV in the case of (a) s u = s d = 0.1 and (b) s u = 0.1, s d = 0. The pink 
short line denotes the SM prediction with the error-bar of 90%C.L. 
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Figure 3: The predicted s u dependence of (a) BR(Kl — > 7r°z/h) and (b) BR(K + —> 7r + z/h) 
at the SUSY scale of 10 TeV, where s d is scanned in the region of 0 ~ 0.3 independent of s u . 
The red dashed lines denote the la experimental bounds for BR(K + —)■ n + vv). The black 
line corresponds to the Grossman-Nir bound together with the experimental upper bound of 
BR(K + —> tt + uu) with 3a. 

and 4 L = sfjf = s d = 0.1 {i = 1, 2) for the up-quark and the down-quark sectors, respectively. 
The Z penguin mediated chargino dominates the SUSY contribution to these branching ratios. 

The SUSY contributions can enhance the branching ratio of Kl —> in eight times 

compared with the SM predictions in Eq.Q, 1.8 x 10~ 10 although it is much smaller than 
the Grossman-Nir bound. On the other hand, the predicted BR(K + —> ir + vv) increases up 
to three times, 2.1 x 10 -10 . It is also noticed that the predicted region of BR(Kl —> n Q vv) 
is reduced to much smaller than 10 -11 due to the cancellation between the SM and SUSY 
contributions. The BR(K + —$■ 7r + z/h) could be reduced to 1.3 x 1CT 11 . 

We discuss the correlation between ex and BR(Kl — > 7r°uu) in Fig. [2l in which (a) 
s u = s d = 0.1 and (b) s u = 0.1, s d = 0. The transverse axis denotes the SUSY contribution 
in \ex\- If the down-squark mixing s d is comparable to the up-squark mixing s u , there is 
no correlation between them as seen in Fig. [2](a), where the Z penguin mediated chargino 
dominates the SUSY contribution of Kl —> n 0 vu, and the gluino box diagram dominates the 
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Figure 4: The predicted BR(K L —> tt°uu) versus BR(K + —y n + uu) at the SUSY scale of 50 
TeV with the mixing angle of s u = s d = 0.3. The pink cross denotes the SM predictions. The 
red dashed lines are the la experimental values for BR(K + —> 7r + uu). The green slanting 
line shows the Grossman-Nir bound [ 13] . 


SUSY contribution of 6k- The gluino contribution of 30% is possible in 6k- 

On the other hand, if the down-squark mixing s d is tiny compared with the up-squark 
mixing s u , the Z penguin mediated chargino dominates both SUSY contributions of A% —y 
7i°vv and 6k- Then, it is found a correlation between them as seen in Fig. [2(b), where the 
chargino contribution to 6k is at most 3%. This correlation is due to the difference of the 
phase structure between the penguin diagram and the box diagram of the chargino. 

In conclusion, 6k could be deviated from the SM prediction in 0(10%) due to the gluino 
box diagram, whereas the Z penguin mediated chargino could enhance the branching ratio 
of Kl —> 7 t°uu from the SM prediction. 

Next, in order to see the mixing angle s u dependence of the branching ratios, we plot the 
predicted regions on BR(Kl — > n°uD) vs. s u and BR(K + —y n + uu) vs. s u planes taking 
s u = 0 ~ 0.3 in FigJH] (a) and (b). We scan s d in the region of 0 ~ 0.3 independent of s u 
although the gluino contribution is much suppressed compared with the chargino one. In 
this plot, the SUSY contribution to 6k is free (0 — 40%), but the experimental constraint of 
|ck-| with the error-bar of 90%C.L. is taken account. We show the upper bound given by the 
Grossman-Nir bound together with the experimental upper bound of BR(K + —y i\ + vv) with 
3cr by the black line, at which the predicted BR(K L —» -k°vv) should be cut. Namely, the 
observed upper bound of BR(K + —$■ ir+vi') gives the constraint for the predicted BR(K L —* 
7 t°uu) at s u larger than 0.2. The precise experimental measurement of BR(K + —> n + vv) will 
lower the predicted upper bound of BR(Kl —> ir 0 vv). 

Let us discuss the case of the SUSY scale of 50 TeV. Fig. [4] shows the predictions on 
the BR(Kl —> and BR(K + —>■ 7 x + vv) plane at the SUSY scale of 50 TeV, where the 

mixing angle is fixed at s u = s d = 0.3. Although the predicted region is reduced considerably 
comparing to the case of the 10 TeV scale in Fig. [H the predicted branching ratio of —>■ 

is enhanced in two times from the SM prediction, and the branching ratio of K + —> tt+vv 
could be enhanced from the SM prediction in three times. 
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Figure 5: The predicted BR(Kl —> tt°uP) versus the SUSY contribution ratio of ex at the 
SUSY scale of 50 TeV in the case of (a) s u = s d = 0.3 and (b) s u = 0.3, s d = 0. The pink 
short line denotes the SM prediction with the error-bar of 90%C.L. 
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Figure 6: The predicted s u dependence of (a) BR(K L —$■ 7 t g vv) and (b) BR(K + —> i t + vv) 
at the SUSY scale of 50 TeV, where s d is scanned in the region of 0 ~ 0.3 independent of s u . 


To see the correlation between ex and the predicted Kl —> n°uP branching ratio, we show 
the branching ratio of Kl —>■ 7i°uu versus the SUSY contribution of ex in Fig. 0 in which 
(a) s u = s d = 0.3 and (b) s u = 0.3, s d = 0. We do not find any correlation between them 
in the Fig. 0(a), where the gluino contribution to ex is still possible up to 10%. However, 
it is found a correlation between them as seen in Fig. 0(b), where the Z penguin mediated 
chargino dominates both SUSY contributions of K L —> n°vv and ex since the down-squark 
mixing Sd vanishes with keeping s u = 0.3. The chargino contribution to ex is at most 2%. 
This correlation is understandable from the difference of the phase structure between the 
penguin diagram and the box diagram of the chargino. 

Thus, even if the SUSY scale is 50 TeV, ex could be deviated from the SM prediction in 
0(10%) due to the gluino box diagram, whereas the chargino process deviates the branching 
ratio of Kl —>■ n°vv from the SM prediction in the factor two. 

Fig0 shows the s u dependence of BR(Kl —> n°uu) and BR(K + —>• n + uu) taking s u = 
0 ~ 0.5 in Fig0 (a) and (b). We also scan s d in the region of 0 ~ 0.3 independent of s u . In 
this plot, the SUSY contribution to ex is free (0 — 40%), but the experimental constraint of 
ex with the error-bar of 90%C.L. is taken account. The predicted BR(K L —> tt°uP) could be 
large up to 8 x 10 —11 , and BR(K + —>■ n + uu) is up to 1.5 x 10~ 10 . Thus, the enhancement 
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Figure 7: The predicted BR(K L —* 7r°e + e _ ) versus BR(K L —>■ n 0 uu) with s u = 0 ~ 0.3 and 
s d = 0 ~ 0.3 at the SUSY scale of 10 TeV. The red solid line denotes the upper-bound of the 
branching ratio BR(Kl —> n 0 e + e~). 


from the SM prediction could be detectable even if the SUSY scale is 50TeV. 

Before closing our numerical study, we would like to discuss correlations to other quantities 
which are sensitive to the NP. They are K L —> 7r°e + e” process and the neutron electric 
dipole moment d n . The K Tj —>- 7T°e + e _ process is induced in similar way to K Tj —>■ 7r°z/P. 
The distinguish feature of Kl —>■ vr°e + e _ mode is the contribution of the photon penguin. 
Moreover, one cannot neglect the long-distance effect from the photon exchange process 
j48| . Thus, the decay amplitude of K L —>■ n°e + e~ has both the short-distance effect and 
the long-distance effect, and the SM prediction of the branching ratio is around 3 x 10 —11 , 
which is comparable to the SM prediction of K L —>• Since our interest here is to check 

whether the SUSY effect does not exceed the experimental bound of K L —y 7r°e + e“, we only 
consider the short distance contribution in our analysis. The experimental bound of the 
branching ratio K L —> 7r°e + e~ is BR(K L —>■ 7r°e + e _ ) exp < 2.8 x 1CT 10 @3J. In the FigH the 
predicted BR(K L 7r°e + e _ ) vs. BR(K L plane are plotted with s u = 0 ~ 0.3 

and s d = 0 ~ 0.3 at the lOTeV scale of the SUSY. There are two predicted lines in this 
figure. Because the decay amplitude A(Kl —> ir°e + e~) is described by the sum of the SM 
and the SUSY contributions, there are two ways of taking the relative phase of ± such 
as A(K l —> 7r°e + e _ ) = A(K L —>■ 7r°e + e^ : SM) ± A(K L —>■ 7r°e + e _ : SUSY), which has 
two solutions giving the same absolute value of the decay amplitude. Then, we have two 
predicted values of BR(Kl — y n 0 e + e~ ) for the certain BR(Ki —y tAvv). The both decay 
processes are dominated by the Z penguin mediated charginos, then, the branching ratios 
are determined by the final state couplings of Zvv and Ze + e~, that is, the weak charges 
Q^zl an d Qzl- Moreover, three flavors of neutrinos are summed for K rj —> Therefore, 

BR(Kl —y 7 t°uu) is significantly larger than BR(Kl —> 7r°e + e _ ). On the other hand, in the 
SM, there are some contributions to Kl —>■ 7T°e + e _ such as the photon exchange processes. 
So, BR(Kl —>■ 7r°e + e _ ) is comparable to BR(Kl vr °vv) in the SM. In conclusion, the 
experimental upper bound of BR(K L — > 7i°e + e ~) excludes the region larger than BR(K L —* 
7 t°uu) = 1.7 x 10” 9 . However, if the long-distance effect is properly included f4S], this 
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constraint becomes somewhat tight or loose depending on the relative sign between the 
SUSY contribution and the long-distance one. 

The neutron electric dipole moment (EDM) d n is well known as the sensitive probe for 
the NP, and so we have studied the correlation between the neutron EDM and the K —y i t°uD 
branching ratio. It is found that our predicted K —» n°uP does not correlate with d n . Suppose 
the SUSY contribution to the chromo-EDM of quarks through the gluon penguin mediated 
gluino [49] -[53], where the left-right mixing term of the down-squark is dominant. In our 
SUSY mass spectra, the left-right mixing is suppressed as discussed in section 3. Moreover, 
the CP violating phase dependence of d n comes from the down-squark mixing matrix whereas 
the phase of K —> n comes from the up-squark mixing matrix. Namely, those phase 
dependences are completely different each other. Therefore, we do not take account of the 
constraint from the experimental upper bound of the neutron EDM in our analyses. 

5 Summary 

We have studied the contribution of the high-scale SUSY to the Kj j —>■ n and K + —y 
n + vv processes by correlating with the CP violating parameter 6k- These rare decays have 
important role of the decision of the CP phase in the CKM matrix, furthermore, they are 
also sensitive to the flavor structure of the NP. 

Taking account of the recent LHC results for the Higgs discovery and the SUSY searches, 
we consider the hight-scale SUSY at the 10 — 50TeV scale. Then, we have discussed the 
SUSY effects to K + —> tt + uD, Kl —> 7r 0 uP and 6k in the framework of the mass eigenstate 
basis of the SUSY particles assuming the non-minimal squark (slepton) flavor mixing. 

We have calculated the SUSY contribution to the branching ratios of Kl —y n and 
K + —> 7where phase parameters are constrained by the observed 6k- The Z penguin 
mediated chargino dominates the SUSY contribution for these decays. At the 10 TeV scale 
of the SUSY, its contribution can enhance the branching ratio of Kl —> n 0 uu in eight times 
compared with the SM predictions whereas the predicted branching ratio BR(K + —y ti + vu) 
increases up to three times of the SM prediction in the case of the up-squark mixing s u = 0.1. 

We have investigated the correlation between 6k and the Kl —>■ tt 0 uu branching ratio. 
Since the gluino box diagram dominates the SUSY contribution of 6k up to 30%, there is no 
correlation between them. However, if the down-squark mixing is neglected compared with 
the up-squark mixing, the chargino process dominates both SUSY contributions of Kl —$■ 
n and 6k- Then, it is found a correlation between them, but the chargino contribution 
to 6k is at most 3%. It is concluded that 6k could be deviated significantly from the SM 
prediction in 0(10%) due to the gluino box process, whereas the chargino process could 
enhance the branching ratio of Kl —> n 0 uu in several times from the SM prediction. 

Our predicted branching ratios depend on the mixing angle s u significantly. The observed 
upper bound of BR(K + —y n + uu) gives the constraint for the predicted BR(Kl —> tt q vu) at 
s u larger than 0.2. 

Even if the SUSY scale is 50 TeV, the chargino process still enhances the branching ratio 
of Kl —>■ 7r°z/P from the SM prediction in the factor two, and the 6k is deviated from the SM 
prediction in 0(10%) unless the down-squark mixing s d is suppressed. 
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We also discuss correlations to the Ki —$■ 7T°e + e~ process and the neutron electric dipole 
moment d n which are sensitive to the NP. 

We expect the measurement of these processes will be improved by the J-PARC KOTO 
experiment and CERN NA62 experiment in the near future. 
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Appendix A : SUSY Spectrum 

In the framework of the MSSM, one obtains the SUSY particle spectrum which is consistent 
with the observed Higgs mass. The numerical analyses have been given in Refs. [MII55]. At 
the SUSY breaking scale A, the quadratic terms in the MSSM potential is given as 


U 2 = m\\Hi\ 2 + m 2 2 \H 2 \ 2 + m 2 (iA • H 2 + h.c .) . 
The mass eigenvalues at the II\ and H 2 = eH 2 system are given 


(34) 


9 772? + mn 

m 2 = —--- =f 

o 


m J — m. 


+ m\ . 


(35) 


Suppose that the MSSM matches with the SM at the SUSY mass scale Q 0 = m 0 . Then, the 
smaller one m 2 _ is identified to be the mass squared of the SM Higgs H with the tachyonic 
mass. The larger one is the mass squared of the orthogonal combination "H, which is 
decoupled from the SM at Q o, that is, m ~ Q 0 . Therefore, we have 


m 2 _ = —m 2 (Q 0 ) , m 2 + = m^(Qo) = ml + ml + m 2 , 


with 


4 / 2 i 2\/ 2 i 2\ 

m 3 = ( m l + m ) ( m 2 + m ) , 

which leads to the mixing angle between Hi and H 2 , (3 as follows: 


(36) 


(37) 


-I— TTX? ~ ~ 

tan 2 (3 = —5 -- , H = cos (3Hi + sin (3H 2 , % — — sin /3Hi + cos (3H 2 ■ (38) 

mg + 

Thus, the Higgs mass parameter m 2 is expressed in terms of m 2 , and tan (3: 

,2 _ 


2 mr — mi tan 2 f3 

772" = 


(39) 


tan 2 (3 — 1 

Below the Qq scale, in which the SM emerges, the scalar potential is the SM one as follows: 


Vsm — ~m 2 \H\ 2 + 


(40) 


Here, the Higgs coupling A is given in terms of the SUSY parameters at the leading order as 

A(Qo ) - AQo) c °t (3 


A(Q„) = hj 2 + 9 ' 2 ) cos 2 2/3 + Hx, 2 (l - Y 


87T 2 


12 


Y, = 


Qo 


, (41) 


and ht is the top Yukawa coupling of the SM. The parameters m 2 and A run with the SM 
Renormalization Group Equation down to the electroweak scale Qew = ^h, and then give 


m 2 H = 2?72 2 (772 j H-) = A(7T2#)u 2 


(42) 
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It is easily seen that the VEV of Higgs, ( H) is v, and (H) = 0, taking account of (Hi) = v cos /3 
and (H 2 ) = v sin /3, where v = 246GeV. 

Let us £x m E = 126GeV, which gives A(Qo) and m 2 (Q 0 ). This experimental input 
constrains the SUSY mass spectrum of the MSSM. We consider the some universal soft 
breaking parameters at the SUSY breaking scale A as follows: 

m Qi( A ) = m u°( A ) = m D^ A ) = A ) = "^( A ) = rn 2 0 (i = 1, 2, 3) , 

Mi (A) = M 2 (A) = M 3 (A) = mi /2 , m\( A) = m 2 (A) = m 2 0 , 

Au(A) = A 0 yu(A) , A D ( A) = A 0 y D ( A) , A E ( A) = A 0 y E ( A) . (43) 

Therefore, there is no flavor mixing at A in the MSSM. However, in order to consider the 
non-minimal flavor mixing framework, we allow the off diagonal components of the squark 
mass matrices at the 10% level, which leads to the flavor mixing of order 0.1. We take these 
flavor mixing angles as free parameters at low energies. 

Now, we have the SUSY five parameters, A, tan/3, m 0 , mi/ 2 , where Q 0 = m 0 . In 
addition to these parameters, we take /i = Qo- By fixing A, Q 0 and tan/3, we tune mi / 2 and 
Aq in order to obtain m 2 (Qo) and Xh(Qo) which realize the correct electroweak vacuum with 
run = 126GeV. Then, we obtain the SUSY particle spectrum. We consider the two case of 
Q 0 = 10 TeV and 50 TeV. The input parameter set and the obtained SUSY mass spectra at 
Q 0 are summarized in Table 1, where we use m t (m t ) = 163.5 ± 2 GeV [43| [56lj . 


Input at A and Q o 

Output at Qo 

at A = 10 17 GeV, 
m 0 = 10 TeV, 
mi /2 = 6.2 TeV, 

A 0 = 25.803 TeV; 
at Q 0 = 10 TeV, 
p = 10 TeV, 
tan/3 = 10 

rrig = 12.8 TeV, m^ = 5.2 TeV, m§ = 2.9 TeV 

rrii L = mi L = 12.2 TeV 

m~ bR = 14.1 TeV, m^ R =8.4 TeV 

n h L ,d L = m c L ,u L = 15.1 TeV 

m s R ,d R — m ~CR,u R = 14.6 TeV, m n = 13.7 TeV 

rrif L = iTh/ TL = 10.4 TeV, rrif R = 9.3 TeV 

m ^ L ,e L = m ^ L ,u eL = 1 °- 8 TeV ; m m,e R = 10-3 TeV 

X t = -0.22, \ H = 0.126 

ml = 1.84857 x 10 8 GeV 2 , m 2 = 1.83996 x 10 6 GeV 2 , m 2 = 8691GeV 2 

at A = 10 lts GeV, 
m 0 = 50 TeV, 
mi /2 = 63.5 TeV, 

A 0 = 109.993 TeV; 
at Qo = 50 TeV, 
p = 50 TeV, 
tan (3 = 4 

mg = 115.6 TeV, m^ = 55.4 TeV, m^ = 33.45 TeV 

m~ bL = mi L = 100.9 TeV 

m bR = 104.0 TeV, m,f R = 83.2 TeV 

m s L J L = m c L ,u L = H0.7 TeV, rn~ R ( j R = 110.7 TeV 

m c R ,u R — 105.0 TeV, m^i = 83.1 TeV 

mf L = m u - TL = 63.6 TeV, mf R = 54.6 TeV 

m idL,e L = m^ L ,v eL = 63.8 TeV, rri jjiR ,e R = 55.0 TeV 

X t = -0.65, A H = 0.1007 

m 2 = 6.49990 x 10 9 GeV 2 , m 2 = 4.06235 x 10 8 GeV 2 , m 2 = 8840GeV 2 


Table 1: Input parameters at A and the obtained SUSY spectra at Q o = 10 and 50TeV. 
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Appendix B : Z penguin amplitude mediated charginos 

We present the expression for the Z penguin amplitude mediated the chargino, -P|l(a ) i n 
our basis m as follows: 


f&txb = Sr E ( r Ci , )h(t'S)J'b/(t4.)J(£ / +)i [>og<° + /2(^,^)] (44) 

J W n a T j 


MUl / r (d)- 


where 


and 


a,/3.I,J 


- 2Sj(U i _)l(U-)';^x , a x' ll f l (x‘ a ,x‘,)-Sp j. , 

(rh« 




= (r^VcKM)? ([/+)? + -^ { RfuV C KM)W 

92 


a 

+ )2 > 




(45) 

(46) 

(47) 


with q = s, d, I = 1 — 6 for up-squarks, and a = 1, 2 for charginos. The Vckm is the CKM 
matrix, and U± are the 2x2 unitary matrices which diagonalize the chargino mass matrix 

M c -. 


U ] _M C U + = -diagMg , (a = 1, 2) . 


(48) 


The fu denotes the yukawa coupling defined by fjjv sin/3 = diag (m u ,m c ,mt). The loop 
integral functions are given as: 


with 


fn(x,y) 


1 


x-y 


x n \ogx 

x — 1 


y n iogy \ 
y-1) ’ 


X 


I 

a 


m: 


m. 


X I 


~ 2 
rrij 

2 ~ ’ 
ho 


where /Uq = Qo is taken in our framework. 


(49) 


(50) 
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